N95- 32432 


IMPLEMENTATION OF A SINGLE-STAGE-TO-ORBIT (SSTO) 
MODEL FOR STABILITY AND CONTROL ANALYSIS 

Final Report 

NASA/ASEE Summer Faculty Fellowship Program— 1994 
Johnson Space Center 


Prepared By: 

Stephen A. Ingalls, MSAE 

Academic Rank: 

Assistant Professor 

College and Department: 

United States Military Academy 
Department of Civil and 

Mechanical Engineering 
West Point, New York 10996 

NASA/JSC 


Directorate: 

Engineering 

Division: 

Navigation, Control, and Aeronautics 

Branch: 

Control and Guidance 

JSC Colleague: 

G. Gene McSwain 

Date Submitted: 

August 5, 1994 

Contract Number: 

NGT-44-005-803 


15-1 



The tendency to combine the three metrics into one will only obscure 
information and should be avoided. The recommendation that each of the 
three be tracked separately is yet another reason to restricting ourselves to 
three metrics. 

It is fully expected that each IPT will avail themselves of sophisticated 
methods to establish the value of the parameters. It may well even be true that 
each IPT uses a different method. However, over time and with experience, it 
is reasonable to suppose that some methods will prevail and become 
common. Here, the intent is to let use determine the method as opposed to 
method determining the use. 

Step 4: Reconciliation 

The top ten plus one list should be compared to the three metrics 
generated to see if the information plays well together. Logic and consistency 
are important. In essence, the IPT has verified their feelings on risk by working 
the problem two different ways, one in paragraph form and the other with the 
metrics. This part of the report should be in paragraph form and discuss how 
the two different parts support each other. 

Step 5: Transition 

Risk management is dynamic. The RM plan requires a periodic review 
to insure it is meeting its required objectives and to insure that the objectives 
themselves have not changed. As the program becomes more operational in 
nature, significant changes in the definition of success will occur forcing 
changes in the IRM plan. Step 5 in the report is a paragraph discussing how 
the office plans to deal with the change from a risk perspective. 

CONCLUSIONS 

As was stated earlier in the paper, this is a tentative beginning of the 
development of the foundations of risk management. Surely significant 
changes and modifications will occur with time. The question now becomes 
one of whether the implementation of such a program is worth the overhead 
that it brings. Can a program afford to implement an IRM plan? 

One of the problems with technical decision making is that non-technical 
people understand neither the principles involved nor the process. Given this, 
they have to trust the technocrat to do what is right and correct One thing is 
resoundingly clear, if the community as a whole loses their faith in the 
technocrats to make the "good" decision, they can shut the whole process 
down. If Congress loses its faith in NASA, they can stop the space station. If 
the public loses its faith in a company, they can force financial ruin. If the 
auditors and accountants lose their faith in the design team they can stop a 
project. 

Given this and given the high risk of many technical decisions, the 
question more properly becomes one of can a company or program afford not 
to institute a formal risk management plan? 
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high-performance aircraft are purposely designed statically unstable, 
increasing their maneuverability. Instabilities are easily controlled using 
feedback: pilot-in-the-loop, automatic flight control systems, or both. 

Dynamic stability is the time-response of an aircraft to a perturbation. 
Dynamic instability is also acceptable in an aircraft, presuming the divergence 
from an equilibrium condition is slow enough. These instabilities are also 
dealt with using human or automatic feedback. 

Moment Trim 

Separate from, yet related to, stability is moment trim. For an aircraft 
to achieve an equilibrium flight condition, moments about the center of 
gravity, both longitudinally and laterally, must balance. This balance is 
dependent on static stability derivatives and control effectiveness, measured 
in terms of control derivatives, useful aerodynamic angles of attack, and 
maximum allowed deflections. An altitude and airspeed where moment trim 
is not achievable is an infeasible flight condition, regardless of the level of 
static or dynamic stability. 
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DISCUSSION 


Static Stability 

Rather than look at every force and moment response to every input 
motion variable, an arbitrary set of static stability criterion is established: 

• Velocity disturbances are initially opposed only by forces. 

• Rotational velocity disturbances are initially opposed only by moments. 

• Angle of attack and sideslip disturbances obtained by interpreting the 
velocity disturbances v and w as (3 = v/U! and a = w/U x are initially 
opposed only by moments. 

Static stability is indicated when the initial response of the vehicle is a return 
toward the equilibrium point. 


The change in side-force coefficient, C Y , with change in sideslip angle, 
P, serves to illustrate the shorthand notation used with stability derivatives: 

3Cy _ „ 

* ' " (.) 

The criteria for static stability are outlined in Table 1. 


TABLE 1.- AIRPLANE STATIC STABILITY CRITERIA 


Forces and 
Moments 

Drag 

+ X-Thrust 


Sideforce 
+ Y-Thrust 


Lift 

+ Z-Thrust 
Roll Moment 
+ Roll-Thrust 


Pitch Moment 
+ Pitch-Thrust 


Yaw Moment 
+ Yaw-Thrust 


Perturbed Variable 
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As current work concentrates on the descent-to-land mission phase (the 
aircraft is in unpowered, gliding flight), thrust forces of Table 1 are not 
considered. The only derivative contrary to our rules is the pitching moment 
response to a linear velocity, C mu . Inclusion is a function of its physical 
significance. 

Static Stability Results 


Where possible, stability derivatives were determined directly from the 
LaRC aerodynamic database by applying a forward differencing method 
between available data points. Curve-fits were applied to data which was 
highly non-linear with the perturbed variable of interest (drag with angle of 
attack and longitudinal forces and moments with u-velocity (Mach number)). 
Results of this analysis are presented in Table 2. 



Note 2: Dihedral effect (Cjp) is unstable (positive) at Mach 0.30 and angles 
of attack less than 12.0 degrees. For Mach = 0.60, the vehicle is statically 
unstable at less than 10.0 degrees a. 


Note 3: WB001 is statically unstable for all flight conditions studied. 
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Moment Trim 


The six degrees of freedom available to an aircraft are represented by 
six equations of motion (conservation of linear and angular momentum). 
Modeling aerodynamic forces as a first-order Taylor Series expansion and 
writing the equations for rectilinear flight yields: 

mgsmy = ~(C Da +C D a + C Da>/ S bf +C D J t +C D J tf )qS 
-/Mg sin 0 cosy = (C r# + C Y fi + + C Taf 8 f )qS 

mg cos 0 cosy = (C K +C L a + C^8 bf +C L& 8 t +C^8 ir )qS 
0 = C lf p+C,J al + C lv 8 f 
° = c «. +c ^ a+ ^J bf -h^S, +C mv 8, 

O^C-'P+C^+C^S, 

( 2 ) 

Further modification to these equations is possible by assuming zero 
bank angle, <J> = 0. This substitution decouples the longitudinal set of 
equations from the lateral and is representative of actual mission profiles 
flown from 10,000 feet to the runway. 

Longitudinal Equations of Motion 
mgsmy = -(C D ,+C 0 a + C Dw S v +C D J,+C Dv S f )qS 
-mgcosy = (C t + C^a + C^S^ +C Lh S^C L J t/ )qS 

0 = C„ +CLa + C„ 5„ + CL 5. + C. 5, 

m o m a m S>f tf € m&j tf 

In these equations, flight path angle (y), angle of attack (a), body-flap 
deflection (5^), elevator deflection (5 e ), and tip-fin deflection (5^) are 
unknowns. Reducing the five unknowns to three is possible by first setting 
the body-flap deflection to zero, typically done for shuttle in the final stages 
of descent to ensure ground clearance at touchdown. Eliminating angle of 
attack dependency is possible by solving the set with LaRC aerodynamic 
characteristics for WB001 at each angle of attack from -10.0 to 20.0 degrees 
(subsonic range of useful a). The solution set: y, 8 e , and 5 tf , represents a 
longitudinal trim condition. 
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Lateral-Directional Equations of Motion 
-mg cosy = (Q, p + C Yut 8 ail + C Yv 8 f )qS 
® = Cl,P + Clw Sail + Cl v S f 

°=C n ,P + C n J aiI+ C n J f (4) 

In these equations, flight path angle (y), sideslip angle (|3), aileron deflection 
(8^, and tip-fin deflection (8^) are not known. Either the longitudinal flight 
path angle solution, or another flight path angle of interest, is substituted here, 
leaving three unknowns. It is important to differentiate tip-fin use in the 
longitudinal and lateral-directional cases. Tip-fins are employed for speed 
brake effect longitudinally and rudder effect laterally. 

A trim solution does not guarantee feasibility. Maximum control 
deflections possible with the vehicle must exceed the solution obtained from 
these equation sets. Further, aerodynamic angles obtained must not exceed 
stall values in either a or (3. 

Longitudinal Trim Results 

A note concerning trim solutions for elevator deflection is warranted. 
WBOOl's control configuration is depicted in the schematic below (Figure 1). 
For these calculations, it was assumed both inboard and outboard panels were 
deflected to trim pitching moment. Obvious advantages in this are smaller 
control deflections (more surface area moving to counter the moment), 
resulting in decreased drag due to control movement and a greater percentage 
of available deflection to address dynamic stability concerns. Control 
derivatives for the elevator in equation sets 2 and 3 above are, therefore, the 
addition of drag, lift, and pitching moment derivatives with respect to both 
aileron and elevon movement. Longitudinal trim results are shown in Figure 
2 . 


Solutions at Mach 0.30 yield trim flight path angles on the order of 12 - 
13 degrees (orbiter is * 20.0°) and angles of attack varying from around 17.8 
degrees at 10,000 feet to 12.6 degrees at sea-level. These calculations 
presume a standard atmosphere. Trim angles of attack and elevator 
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deflection trend smaller at higher Mach number, while flight path angle 
increases to values in excess of 20 degrees at 10,000 feet. These data bound 
possible trim solutions and this effort does not design a final descent 
trajectory. 



FIGURE 1.- WB001 CONTROL 
ORGANIZATION SCHEMATIC 


Lateral-Directional Trim Results 

Lateral-directional trim was 
attempted in a similar manner. Rather 
than look at conditions throughout the 
final descent, trim was attempted for a 
design cross-wind condition at landing 
of 25 knots. Assuming touchdown 
velocities comparable to shuttle (= 
200 knots), this results in a design 
sideslip angle of approximately 7 
degrees. 


Longitudinal Trim Solutions 

WB001 @ xcg = 0.714 
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FIGURE 2.- LONGITUDINAL TRIM SOLUTIONS (a, y, and 8 e ) 
FOR TWO MACH NUMBERS VS. ALTITUDE 
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Further, orbiter design vertical velocities at touchdown are 
approximately 2.0 feet/second. Using this and the assumed touchdown 
velocity of around 200 knots results in a flight path angle of approximately 
zero (0.34°) for substitution into the lateral equation set. 

Attempts to obtain a trim solution to the lateral equations of motion 
proved impossible, indicating the vehicle incapable of landing in the design 
cross-wind. This, coupled with C n p static instabilities, suggested an analysis 
to improve the SSTO's inherent lateral stability. 


Tip-Fin Sizing 

Early landing simulation of the SSTO indicated lateral divergence from 
design trajectories in both roll and yaw with small (3 (Figure 3). Further 
attempts to land the vehicle indicated increasing lateral force and moment 
effects resulted in more favorable performance. Discussions with engineers 
at LaRC indicated a pending effort to reevaluate tip-fin size. 


Sizing lateral control surfaces to produce a desired C np derivative value 
is common practice in the preliminary design of aircraft. Considering vehicle 
weathercock stability (C n p) as a superposition of effects from the wing-body 
and lateral control surfaces, in this case, the tip-fins, yields: 


C n =C„ +C„ 

n p_ttp-ju, 


(5) 


The tip-fin contribution is further expressed as: 


C n =C 


da. 


A 30 Hr 


S b 


( 6 ) 


where. 


^Latip-fin 

= Tip-Fin Lift-Curve Slope 


= Change in Sidewash with Sideslip Angle 

ntf 

= Ratio of Dynamic Pressure at the Tip-Fin to the Wing 

s tf 

= Tip-Fin Surface Area [ft 2 ] 

S 

= Wing Area [ft 2 ] 

x tf 

= Distance from Center-of-Gravity to Tip-Fin Aerodynamic 
Center [ft] 

b 

= Wing Span [ft] 
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FIGURE 3.- SHUTTLE ENGINEERING SIMULATOR STRIP PLOTS COMPARING 
ROLL RATE PEG/SEC], PITCH RATE pEG/SEC] AND SIDESLIP ANGLE PEG] 
VS. TIME [SECONDS] FOR THE SSTO TO ORBITERIN STS-64 

CONFIGURATION 

If we neglect sidewash dependence on p (9a/9P) and assume the 
dynamic pressure ratio (r) v ) is 1, equation 6 reduces to: 

S, X,, 

Q _ Q £ £ 

n PJp-M O L 

O D (7) 

C La _tip-fin was calculated as 2.174/radian using methods outlined in the 
United States Air Force Stability and Control Data Compendium 
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(DATCOM). Comparison of LaRC tip-fin yaw moment to tip-fin side force 
at the same deflection resulted in a value for X tf (moment arm of the tip-fin 
control surface to vehicle center-of-gravity) of 56.112 feet. Substituting a 
wing area of 4,211 ft 2 and span equal to 95.72 ft, equation 7 shows how 
increasing tip-fin surface area results in a linear increase in C n p tip-fin- 

A vehicle C n p of 0.100/radian was selected and is representative of the 
orbiter's subsonic lateral stability in this derivative. WBOOl's current tip-fin 
surface area is 214.95 ft 2 . Increasing area to obtain this degree of inherent 
stability requires a growth factor of 6.78 (total lateral control surface area 
required increases to 1,457 ft 2 . 

This offers one method of improving vehicle lateral stability and cross- 
wind landing capability. Enhancing basic wing-body weathercock stability or 
employing alternative methods of improving tip-fin effectiveness without a 
size increase are equally viable options. 
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CONCLUSIONS 


This study concludes the WBOOl's current level of lateral stability is 
insufficient. In particular, lateral trim solutions for desired cross-wind 
components in the landing environment were not possible. Instabilities in roll 
and pitch 'tuck' also exist, but appear less significant than deficiencies in the 
C n p derivative. Longitudinal trim solutions indicate steady-state glide path 
angles on the order of 12 degrees at Mach = 0.30, increasing to slightly 
greater than 20 degrees at M = 0.60. One attempt at improving the vehicle's 
lateral instabilities looked at sizing lateral control surfaces to a prescribed 
degree of inherent stability. A value comparable to subsonic stability of the 
orbiter (C n p = 0.100/radian) was selected and required increasing tip-fin 
surface area 6.78 times (from around 215 to 1,457 ft 2 ). This carries a 
commensurate increase in aircraft weight, highly undesirable in a SSTO 
configuration. Other methods of enhancing the vehicle's lateral stability are 
available and require further study. 



REFERENCES 


1 . Access to Space. NASA Advanced Technology Team (Option 3), 

Volume 3. Final Report. 1993. 

2. Engelund, W. "SSV-001, 6-DOF Aero Database." Preliminary 

Aerodynamic Database. March 1994. 

3. Hoak, D. E. et al. "USAF Stability and Control DATCOM." Data 

Compendium, 1978. 

4. Lan, C. E. and J. Roskam. Airplane Aerodynamics and Performance . 

Roskam Aviation and Engineering Corporation, Ottawa, Kansas, 
1988. 

5. Lepsch, R. "Weight Statement - Level III, Unmanned SSV, Duel-Fuel, 

RD-701." Data Table, April 1994. 

6. Rockwell International. "Operational Aerodynamic Data Book, 

Orbiter Vehicle Aerodynamic Data, Volume 3." Rockwell 
International, 1985. 

7. Rockwell International. "Space Shuttle Ascent GN&C Data Book 

Appendix." Contract NAS9-1 8500. September 1992. 

8. Roskam, J. Airplane Flight Dynamics and Automatic Flight Controls. 

Roskam Aviation and Engineering Corporation, Ottawa, Kansas, 
1982. 

9. Springer, A. "Initial Ascent and Entry Aerodynamic Database for the 

Wing/Body 001 (WB001) Reusable Launch Vehicle (RLV) 
Configuration." Memorandum, July 1994. 


15-14 



N95- 32433 


3 V— 

•7 “ f 
/ ' 1 


SCANNERS, OPTICAL CHARACTER READERS, CYRILLIC 
ALPHABET AND RUSSIAN TRANSLATIONS 


Final Report 

NASA/ASEE Summer Faculty Fellowship Program 1994 
Johnson Space Center 


Prepared By: 

Academic Rank: 

College and Department: 

NASA/JSC 

Directorate: 

Division: 

Branch: 

JSC Colleague: 

Date Submitted: 

Contract Number: 


Gordon G. Johnson, Ph.D. 
Professor 

University of Houston 
Department of Mathematics 
Houston, Texas 77204 

Information Systems 
Technical Systems 
Client Server Systems 
Lui Wang 
August 16,1994 
NGT-44-005 -803 


16-1 



